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This letter deals with triply-periodic (crystalline) solutions in a family of Skyrme systems, namely where the field takes 
values in the squashed 3-sphere. The family includes the standard Skyrme model (round 3-sphere), and the Skyrme- 
Faddeev case (maximal squashing). In the round case, the lowest-energy crystal is the well-known cubic lattice of 
half- skyrmions; but in the squashed case the minimal-energy crystal structures turn out to be different. We describe 
some of the solutions that arise, including arrays of vortices and multi-sheeted structures. 
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1. Introduction 

In the basic SU(2) Skyrme model, the solution with the 
lowest energy-per-charge En is the triply-periodic skyrme 
crystal |1 . For a given value of the topological charge 
(winding number, or baryon number) TV, there are in gen- 
eral many isolated A^-skyrmion solutions, i.e. local min- 
ima of the energy; but the belief is that for N ^> 1, the 
lowest-energy solution will resemble a chunk of the skyrme 
crystal. 

However, in the Skyrme-Faddeev system, where the 
field takes values in S 2 and the topological charge N is 
the Hopf number, the situation is less clear. If there were 
a triply-periodic crystalline solution with En = c (where 
c is a constant independent of A/"), as in the Skyrme case, 
then a large crystalline chunk would have En ~ c. But 
we know [2] that there exist solutions with En ~ TV -1 / 4 , 
so such a crystal chunk certainly could not be the global 
energy minimum among fields of charge TV, although it 
could be a local minimum. In general, it seems to be the 
case that Hopf solitons tend to clump into a tangle; various 
examples, for low values of TV, may be seen in [3 , 4 , 5j|6j[7]- 

The aim of this letter is to explore this difference in 
behaviour, by investigating triply-periodic crystalline so- 
lutions in a system which interpolates between the two 
cases above. This generalized Skyrme system was intro- 
duced in [8]; it is labelled by a parameter a G [0, 1], in such 
a way that a = gives the 'pure' Skyrme model, whereas 
a = 1 gives the Skyrme-Faddeev model. In particular, 
we wish to see what happens to the skyrme crystal as the 
system deforms away from the basic Skyrme case. 

The interpolating system is very natural geometrically. 
For simplicity, let us restrict to static fields in all of what 



follows. Let E be a compact 3-dimensional, oriented, con- 
nected Riemannian manifold. Then for fields <£> : R 3 — >• E, 
the Skyrme energy functional E[<&\ has a natural defini- 
tion, depending in particular on the metric of E [9 j. With 
the usual boundary condition $(x J ) — > <3>o as | — » oo, 
where <3>o is some specified point on E, the topological 
charge N is defined to be the degree of the map <£. Our 
system is obtained by taking E to be the squashed 3- 
sphere (Berger sphere): the standard 3-sphere squashed 
along its Hopf fibres by a factor 1 — a. If a = 0, there 
is no squashing: the target space is just the round sphere 
SU(2), and we have the Skyrme model. In the degenerate 
case a = 1, the 3-sphere becomes a 2-sphere, the winding 
number N becomes the Hopf number, and we have the 
Skyrme-Faddeev model. 

The system has a potential condensed-matter inter- 
pretation in which the two complex fields (^1,^2) rep- 
resent two flavours of Cooper pairs [10]. The parameter 
a then appears as the coefficient of a term J M J M , where 
J M = d^Z is the current density. In particular, arrays 
might be of interest in this context. 



2. The Generalized Skyrme System 

In this section, we give some details of the one-parameter 
interpolating family of Skyrme systems. The field is de- 
noted = ($1, <£>2, $3, $4), where each component &p is 
a function of the spatial coordinates x J = (x x ,x 2 ,x 3 ) = 
(x,y,z), and the constraint ^^^^ = 1 is imposed. The 
energy density for the system, parametrized by a G [0, 1], 
is 



£ = X 2 [(d j ^)(d j ^)-aP j P j }+X 4 [2(l-a)F^ F> ^+aQ*&] 

(1) 
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with f2/3 7 being a symplectic form with non-zero compo- 
nents Vt\2 — — ^21 = — ^34 = ^43 = 1- The coupling 
constants A2 and A4 can be scaled as desired, by chang- 
ing the units of energy and length; in what follows, we 
shall use the 'geometrical' choice A2 = l/[47r 2 (3 — a)] and 
A 4 = 1/[4tt 2 (3 - 2a)] as in [8 . The effect of this is that 
the energy of a 1-skyrmion on R 3 is approximately inde- 
pendent of a {E 1 « 1.22). 

If we are thinking of the Skyrme model as a nonlin- 
ear theory of pions, then the deformation when a becomes 
positive may be viewed as putting one of the three pions on 
a different footing from the other two. To see this explic- 
itly, we impose the boundary condition ^ — >> (0,0,0,1) 
as | — > 00, and identify ($1, $2, $3) as the pion fields in 
the usual way. Then in the asymptotic region the energy 
density ([!]) becomes 

£ « \ 2 [(d j $ 1 ) 2 + (d^ 2 ) 2 + (1 - a)(d^ 3 n (2) 

If a = 1, then one of the fields becomes non-dynamic. 
This amounts to passing from the original target space 
S 3 ^ SU(2) to the quotient space S 2 SU(2)/U(1), by 
the Hopf projection 



^1 = 2($i$ 3 - $ 2 $ 4 ) 
^2 = 2(^2^3 + ^1^4) 

^3 = 



■*2 



(3) 



The relevant field is then the unit 3- vector ip, and the 
energy density Q becomes 



£ 



32tt 2 



(^) 2 + -(G ifc ) 2 



(4) 



where Gjk = ^-(djip) x (d^)- This is the Skyrme-Faddeev 
energy [TT] . 

The energy is the integral of |l]) over R 3 , or 

over a fundamental cell if one is dealing with periodic 
fields. There is a Bogomolny-type lower bound [9] on 
the energy of fields with topological charge iV, namely 
E > 6A/VA2A4 vol (E). Using the fact that the volume of 
the squashed 3-sphere is 27r 2 Vl — a, and with our choice 
of A2 and A4, this gives 



E N > 



y/1 — a 



-a/3)(l -2a/3) 



(5) 



In the limit a = 1, the bound ([5| becomes trivial. 
There are various nontrivial bounds in this case. For iso- 
lated Hopf solitons on R 3 , we know that En > c7V _1//4 , 
where c is a constant; it is conjectured that this bound 



holds for c = 1 [12]. This bound is not valid for triply- 
periodic fields : T 3 — >> £ 2 ; but in that case there is 
another lower bound, which follows from the general for- 
mula given in [13J: if is triply-periodic, and L is the 
largest of its periods, then its normalized energy satisfies 
En > tt/(2L). 

For a = 0, the expression ([!]) gives the standard static 
Skyrme energy, and it has an obvious 0(4) symmetry. For 
a ^ 0, the symmetry is broken to U(2), this being the 
subgroup of 0(4) which preserves the symplectic form ft. 

In this letter, we are focusing on crystalline configura- 
tions, so the fields are periodic in (x 1 ,^ 2 ,^ 3 ) with periods 
(^1,^2,^3) respectively. In effect, $ is a map from T 3 to 
5 3 , and it is classified topologically by its degree N G Z. 
The topological classification of maps from T 3 to S 2 , which 
is relevant to the a = 1 case, is more complicated. But 
if we restrict to fields which are algebraically inessential, 
meaning that the 2-form Gjk belongs to the trivial coho- 
mology class in i/ 2 (T 3 ,Z), then the classification remains 
a single integer N G Z (cf. [14 j). The Hopf charge of the 
Hopf projection ip of $ is equal to the degree of <£>, so it is 
consistent to let N denote either. 



3. Crystals in the Generalized Skyrme Family 

In this section, we investigate various triply-periodic 
solutions in the family of generalized Skyrme systems para- 
metrized by a. There are many of these, corresponding to 
different local minima of the energy, and we shall describe 
only some of them. 

The results that follow are obtained by numerical min- 
imization of the energy functional E[$\. More specifically, 
the expression for the energy is replaced by a second- 
order finite-difference approximation on a rectangular lat- 
tice with periodic boundary conditions, and the minima 
are found by conjugate-gradient relaxation. Testing with 
various different values of the lattice spacing indicates that 
the accuracy, for example in the energy E, is better than 
1%. As remarked above, there are many local minima, and 
the initial configuration determines which one of these is 
obtained after flowing down the energy gradient. We also 
vary the side- lengths Li, L2 and L3 of the fundamental 
cell, so as to get the lowest possible value of E[<&\. 

In the original Skyrme system (a = 0), the lowest- 
energy crystal is a cubic lattice [1 . Its unit cell contains 
eight half-skyrmions, and has N = 4. We begin by taking 
this as the initial configuration, and minimizing E[<$>] for 
various values of a. 

In the a = case, the minimum occurs at L\ = Lo = 
L3 = 4.7; but the side- lengths cease to be all equal when 
a > 0. This occurs because the system has less symmetry 
than when a = 0; in particular, one of the spatial direc- 
tions, which we take to be the x 3 -direction, is now on a 
different footing from the other two. For a > 0, the half- 
skyrmions become elongated in the x 3 -direction, and the 
optimal side-lengths satisfy L\ = L 2 > L 3 . These aspects 
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are illustrated in Figure 1. The upper- left plot shows the 
optimal (minimal-energy) values of L\ = L2 and L3, as 
functions of a, in the range < a < 0.9. We see that in 
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Figure 1: The V+AV+V+AV solution. Optimal values of Li = L 2 , 
and L3; and the energy density for a = 0.3, 0.6, 0.9. 

the Skyrme case a = 0, we have L\ = Lo = L3 = 4.7, as 
noted above; but as a increases, the lattice cells are com- 
pressed in the x 3 -direction and expanded in the x 1 - and 
£ 2 -directions. The behaviour in the a = 1 limit will be 
discussed in more detail below. The other three pictures 
in Figure 1 plot the energy density £, or more precisely 
the surfaces where 8 equals 0.7 times its maximum value, 
for three values of a ranging from 0.3 to 0.9. We see that 
as a increases, the half-skyrmions join up pairwise in the 
x 3 -direction, to form four parallel vortices; these are in 
fact two vortices and two antivortices, and so this field 
is called the V+AV+V+AV solution. In the a 1 limit, 
the vortices become homogeneous in the x 3 -direction; such 
vortices will be discussed in more detail below. 

The symmetries of the a > fields depicted in Figure 1 
are a subset of the a = skyrme-crystal symmetries [1 , 
and are generated by 

x j ^ Lj - x\ $4 ^ -$ 4 for j = 1, 2, 3; (6) 

x j ^ Lj/2 - x\ <S>j i-> -<S>j for j = 1, 2, 3; (7) 

(yX , x , x ) 1 y (yX , x , x ) , 

($1,$2,$3,$4) ^ ($ 2 ,$l,$3,-$4). (8) 

The normalized energy of these solutions is plotted in 
Figure 2, as a function of a. The other plots in Figure 2 are 
the Bogomolny bound ([5|, and the energies of two other 
triply-periodic solutions which will be described below. 
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Figure 2: Normalized energy of crystalline solutions, versus a. 

At this point, let us say more about x 3 -homogeneous 
vortices. The field of a p- vortex with unit x 3 -twist, located 
on the x 3 -axis, can be put in the form 

$1 + i$ 2 w (x 1 + ix 2 ) p , $3 + i$ 4 * exp (27rix 3 /L 3 ) (9) 

near that axis. If p > 0, we refer to it as a p-vortex, 
whereas if p < we refer to it as a p-antivortex. If we 
restrict the Hopf projection i/j of § to the plane x 3 = 0, we 
get a map from T 2 to 5 2 ; the degree of this map is the total 
vortex number in the x 3 -direction (taking the multiplicity 
p of each vortex into account). This total vortex number 
has to be zero, since it is the degree of a Hopf projection. In 
other words, there have to be an equal number of vortices 
and antivortices. The topological charge N is the sum of 
the absolute values of the vortex numbers; in other words, 
vortices and antivortices both contribute positively to N. 
The a = 0.9 field illustrated in Figure 1 is close to an N = 
4 multi-vortex configuration, with two 1-vortices and two 
1-antivortices, and hence is referred to as V+AV+V+AV. 

It is worth noting that analogous V+AV+V+AV con- 
figurations can appear in the a = case, for values of the 
side- lengths Lj which differ from the optimal ones [15]. 

Skyrme vortices, or more precisely vortex- antivortex 
pairs, were discussed in [16]. In particular, it was shown 
there that if one has a parallel vortex-antivortex pair sep- 
arated by a (large) distance D, then there is an attractive 
force between them, the leading term of which is propor- 
tional to 1/D. In our a-family, the expression for the 
attractive force is the same, except that it acquires a fac- 
tor of 1 — a. In other words, vortices and antivortices 
attract one another, as long as a < 1. In the a —> 1 limit, 
however, the leading-order term vanishes. It is not known 
whether there is still an attractive vortex-antivortex force 
in this limit, but numerical experiments suggest that the 
force becomes repulsive. In particular, the two vortices 
and two antivortices in our picture repel each other when 
a = 1, and therefore L\^L^ —> 00 as a — » 1, consistent 
with the upper-left plot in Figure 1. 

The second vortex-type solution featuring in Figure 2 
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is 2V+2AV, and it consists of a 2-vortex and a parallel 
2-antivortex. Its energy is plotted for a > 0.6, and it is 
apparent that in this range, its energy is less than that of 
V+AV+V+AV. In this case, the numerical evidence again 
suggests that L±, L 2 —> 00 as a —> 1, implying that the 2- 
vortex and the 2-antivortex repel each other in that limit. 
The four pictures in Figure 3 plot the energy density iso- 
surfaces where £ equals 0.8 times its maximum value, for 
four values of a ranging from 0.6 to 0.99. 
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Figure 3: The 2V+2AV solution: energy densities for a = 
0.6, 0.7, 0.8, 0.99. 

Note that in the a = 1 (Skyrme-Faddeev) system, we 
can have vortices without antivortex partners — the re- 
striction of the net vortex number being zero does not 
apply in this case. Such Hopf- vortex fields have been the 
subject of several studies, such as [HJ[5j[7]. The results 
described above suggest that in the Skyrme-Faddeev sys- 
tem, the normalized energy of a p-vortex decreases as p 
increases. This indeed turns out to be the case, as we now 
explain in more detail. 

Consider fields if with energy density Q on R 2 x S 1 . So 
i/j is periodic in x 3 with period L, and it satisfies the bound- 
ary condition ip — ^ (0,0, 1) as p — » 00, where x 1 + ix 2 = 
pe l6 . More specifically, consider rotationally-symmetric p- 
vortices centred on the x 3 -axis: these will have the form 

ipi + i^2 = sin (/) exp(ip# + 2irix 3 /L), ijj 3 = cos (/), 

(10) 

where / = f(p) satisfies /(0) = tt and f(p) — > as p — > 00. 
Finding the minimal energy E(p) for various values of the 
vortex number p is a straightforward numerical computa- 
tion, and this was done for p = 1, . . . , 5. Let us, as usual, 
normalize the energy E(p) by dividing it by the topologi- 
cal charge p. The result is that E(p) depends linearly on 
1/p: 

E(p) « 0.338/p + 0.551; 

and it is this value (with, respectively, p = 1 and p = 2) 
which is used for the a = 1 ends of the V+AV+V+AV 
and 2V+2AV plots in Figure 2. 



In view of the above, it seems likely that triply-periodic 
solutions with even lower values of En than V+AV+V+AV 
and 2V+2AV could be constructed by assembling parallel 
arrays of p- vortices and p-antivortices for p > 3; but this 
has not been investigated further. 

The solutions corresponding to the final plot in Fig- 
ure 2, referred to as "multi-sheet", have a lower value of 
energy-per-charge than the others mentioned previously if 
a > 0.9, and do not resemble parallel vortex- antivortex 
configurations for a close to 1. Figure 4 shows the energy 
density isosur faces where £ equals 0.8 times its maximum 
value, for four values of a ranging from 0.6 to 0.99. As 
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Figure 4: The " mu It i- sheet" solution: energy densities for a = 
0.6, 0.8, 0.9, 0.99. 

a approaches 1, they are closely approximated by fields 
which are homogeneous in x 1 and x 2 , in fact of the form 

Z 1 := $! + i$ 2 = sin(/) exp(27ri£ 2 x 2 /^2) , (11) 
Z 2 := $4 + i$3 = cos(/) exp(27ri£:ix 1 /^i) , (12) 

where S\ = ±1, 82 = ±1, and / = f(x 3 ). Their energy 
density depends only on x 3 , and is peaked on sheets or- 
thogonal to the x 3 -axis. Such fields arise if we impose the 
symmetries generated by 

x 1 ^ x 1 + c, (Zi, Z 2 ) (Zi, Z 2 exp(27ri£:ic/L 1 )); (13) 

x 2 ^ x 2 + c, (Zi, Z 2 ) 1 y (Zi exp(27ri£ 2 c/L 2 ), Z 2 ); (14) 

(x\x 2 ) ^ (-x\ -x 2 ), (Zi, Z 2 ) ^ (ZT,^). (15) 

Note that such transformations on (Zi,Z 2 ) preserve the 
energy ([!]). In other words, translations in x 1 or x 2 , and 
an inversion, can be compensated by a symmetry trans- 
formation of the field. 

In order to have periodicity in x 3 , and for the topo- 
logical charge N to be nonzero, we need to arrange things 
rather carefully. A periodic field with N = 4 is obtained 
by imposing f(kL 3 /4) = kir/2 for k = 0,1,2,3,4, and 

' (1,1) for < x 3 < L3/4, 

I (-1,1) for L 3 /4<x 3 <L 3 /2, 

l1 ' 2j "I (-1,-1) for L 3 /2<x 3 <3L 3 /4, 

k (1,-1) for 3L 3 /4 < x 3 < L 3 . 



Note that the resulting field is continuous. 

Substituting ([TT][l2| into the formula for the energy 
density gives an expression of the form £ = A(f) (f) 2 + 
B(f), where /' = df/dx 3 . Then a standard Bogomolny 
argument implies that, for given values of L\ and L2 (and 
a), the normalized energy E attains a minimum value 

r7v/2 



E{a,L 1 ,L 2 ) = 2L 1 L 2 f VABdf, (16) 
Jo 

when L3 is given by 



As =4 



JO 



Bdf. 



(17) 



The next step is to minimize ( 16 ) with respect to L\ and 
L2, which was done using numerical integration. The re- 



sulting minimal energy for the ansatz ( [TTj 12) approaches 
that of the multi-sheet solution (which is not quite homo- 
geneous in x 1 and x 2 ) as a — > 1. For example, if a = 0.98, 
then the ansatz energy (16) is only 0.6% higher than that 



of the actual multi-sheet solution, and the fields are almost 
identical. 

For a > 0.9, the multi-sheet solution has the lowest 
energy-per- charge of those that we investigated. However, 
the a —> 1 limit of this solution, and of the ansatz, is rather 
pathological, since the optimal values of L\ = L2 tend to 



zero. In fact, the energy (16) equals 



E = 2ir 



tt/2 , 

y 2A| sin 2 (2 f)L\ + 47r 2 A 2 A 4 sin 4 (2/) df 

(18) 

when a = 1 (and where we have put L\ = L2), from 
which it is clear that its minimal value is attained when 
L\ = 0. At first sight, this behaviour may seem paradox- 
ical, since the usual view is that the fourth-order Skyrme 
term in the energy prevents an object with nontrivial 3- 
dimensional topology from shrinking to zero volume. This 
view is indeed correct for < a < 1: in particular, if a 
lies in this range, then a triply-periodic configuration with 
nonzero topological charge TV has an energy which diverges 
if any of its periods tends to zero. But in the a = 1 limit, 
namely for the Skyrme-Faddeev system, this property no 



longer holds. One can see directly from (11, 12), or rather 



its Hopf projection, why this is so. If one scales each of the 
periods Lj , then each term in the expression of the energy 
of ip scales differently. In particular, the term J (G12) 2 has 
the form 



(Gi 2 ) 2 dx l dx 2 dx 3 = K 



L1L2 ' 



where K is its value when L\ = L2 = L% = 1. This is 
the only term which prevents L\ and L2 from going to 



zero simultaneously. But for (11, 12) the quantity G12 is 



identically zero, and consequently one can always lower 
the energy by reducing L\ and L2. 



4. Conclusions 

We have studied triply-periodic stable solutions (local 
minima of the energy) in a family of Skyrme-type systems 
parametrized by a G [0,1], which interpolates between 
the standard Skyrme model (at a = 0) and the Skyrme- 
Faddeev system (at a = 1). At a = 0, the lowest-energy 
crystal resembles a cubic lattice of half-skyrmion particles, 
and this picture persists near a = 0, as one would expect. 
But for larger values of a, various other configurations 
are preferred (that is, have lower energy-per-charge): for 
example, vortex- antivortex arrays and multi-sheet struc- 
tures. 

All these structures are somewhat problematic in the 
a 1 limit, and it remains unclear whether or not there 
exists a smooth Skyrme-Faddeev crystal which is at least 
stable under small perturbations. As remarked in the in- 
troduction, the sublinear behaviour of the energy perhaps 
makes this unlikely. But the possibility remains open, and 
is worth further investigation. 
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